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ABSTRACT 

Ocean planets are volatile rich planets, not present in our Solar System, which are 
thought to be dominated by deep, global oceans. This results in the formation of 
high-pressure water ice, separating the planetary crust from the liquid ocean and, 
thus, also from the atmosphere. Therefore, instead of a carbonate-silicate cycle like on 
the Earth, the atmospheric carbon dioxide concentration is governed by the capability 
of the ocean to dissolve carbon dioxide (CO2). 

In our study, we focus on the CO2 cycle between the atmosphere and the ocean 
which determines the atmospheric CO2 content. The atmospheric amount of CO2 is 
a fundamental quantity for assessing the potential habitability of the planet’s sur¬ 
face because of its strong greenhouse effect, which determines the planetary surface 
temperature to a large degree. In contrast to the stabilising carbonate-silicate cycle 
regulating the long-term CO2 inventory of the Earth atmosphere, we hnd that the 
CO2 cycle feedback on ocean planets is negative and has strong destabilising effects 
on the planetary climate. By using a chemistry model for oceanic CO2 dissolution 
and an atmospheric model for exoplanets, we show that the CO2 feedback cycle can 
severely limit the extension of the habitable zone for ocean planets. 

Key words: astrobiology, planets and satellites: atmospheres, planets and satellites: 
terrestrial planets, planets and satellites: oceans 


1 INTRODUCTION 


The question of habitability for exoplanets is usually linked 
to the concept of the so called habitable zone (HZ), i.e. 
the range of orbital distances around a host star, where a 
rocky/terrestrial planet can in principle maintain liquid wa¬ 
ter on its surface over an extended period of time. Liquid 
water is - to our current knowledge - one of the most impor¬ 
tant requirement for life as we know it (e.g. Lammer et al. 
2009). However, the classical HZ as introduced by Kasting 


et al. (1993) is strictly dehned for an Earth-like planet, i.e. 


for a planet with a partially rocky surface, a water content 
of one Earth ocean, plate tectonics, and with all the geolog¬ 
ical cycles (in particular, the carbonate-silicate cycle) as on 
Earth. Given the diversity of the already known exoplanets, 
it is highly unlikely that the majority of the low-mass plan- 
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ets are entirely Earth-like and can be described within the 
concept of the classical HZ. 


One particular, highly interesting class of exoplanets 
are ocean planets (or waterworlds), i.e. planets with a much 
higher water content than Earth. Several candidates for mas¬ 
sive water-rich planets have already been found. Due to their 
small orbital distances, none of these can be considered to be 
a potential habitable planet, though. The most well-studied 


example so far is the transiting super-Earth GJ1214 b (Ghar- 
bonneau et al.|[2009 ). Its mean density derived from obser¬ 


vations suggests that roughly half the planet’s bulk compo¬ 
sition could be made of water. In general, the larger radii 
of ocean planets compared to more rocky planets of equal 
mass (Leger et al. 2004 Alibert 2014| , would make them 
very interesting targets for planet detection and character¬ 
ization with, for example, CHEOPS (GHaracterising ExO- 
Planet Satellite, Fortier et al. |2014| , TESS (Transiting Ex- 
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oplanet Survey Satellite, Ricker et al.|2014 |, PLATO (Plan- 
etary Transits and Oscillations of stars, Rauer et al. pMH, 


or the JWST (Janies Webb Space Telescope, Gardner et al. 


2006). Being completely covered by a deep water envelope. 


they fall within the apparent gap between terrestrial plan¬ 
ets with a (partially) rocky surface and the gas giant planets 
with a high amount of hydrogen and helium. 

In this study we focus on the potential habitability of 
low-mass deep ocean planets. First considered by |Leger et al.| 
(2004), these types of planets would form beyond the ice line 


within a protoplanet ary disk, where they collect a large frac¬ 
tion of water as well as other volatiles as ices and later mi¬ 
grate inwards into the habitable zone ( Kuchner|2003 ). These 
considerations are supported by the results from planetary 
formation models, where migrated terrestrial low-mass plan¬ 
ets with a large water inventory (> 100 Earth oceans) lo¬ 
cated in the classically defined habitable zone (HZ) are a 
very common result ( Alibert|2014[ Alibert et al.| 20 f^ Thi- 
abaud et al.||2014 ) 


Studies on the internal structure of such planets by 


Leger et al. (2004) and Alibert (2014) suggest that they 
will have ocean depths in excess of 100 km, which makes 
the presence of a continental landmass highly unlikely. Cat¬ 


egorized as Class IV habitats by hammer et al. (2009) these 


planets can offer one of the most important requirements 
for life: liquid water. Of course, such a large reservoir of 
water (H 2 O) also influences climate processes. One impor¬ 
tant consequence is, for example, the formation of high- 
pressure water ice (ice phases VI and VII) at the bottom 
of the ocean, which prevents the immediate contact of the 
planetary crust with the liquid ocean. Thus, the well-known 


carbonate-silicate cycle (Walker et al. 1981 Kasting et al. 


1993) is unable to operate as a long-term climate stabilisa¬ 
tion mechanism regulating the atmospheric carbon dioxide 
content. Abbot et al. ( 2012 ) also showed that below a land- 


mass fraction of one percent, the usual seafloor weathering, 
which is an important part of the carbonate-silicate cycle, 
becomes inoperative, such that the stabilising weathering- 
climate feedback on ocean planets does not work even if the 
liquid ocean is still in contact with the planetary crust. 

Since the ocean planet originates from beyond the ice 
line in the protoplanetary disk, the total amount of CO 2 ac¬ 
cumulated during its formation is also thought to be much 
higher than the Earths inventory (see e.g. Leger et al. (2004), 


Alibert (2014), or Marboeuf et al.| ( [2014 )). As discussed m 
Leger et al. (2004), the CO 2 amount present above the high- 


pressure water ice layer may, however, be strongly limited by 
gravitational sequestration during the planetary formation 
(i.e. most of the CO 2 is trapped below the water ice layer), 
clathrate formation, or atmospheric loss. In case of negligi¬ 
ble atmospheric loss processes of CO 2 or very inefficient for¬ 
mation of CO 2 -H 2 O clathrates, the total (initial) inventory 
of CO 2 present in the ocean and the atmosphere remains 
constant over time, because CO 2 cannot interact with the 
planetary interior via subduction or outgassing. Then, the 
partial pressure of CO 2 in the atmosphere is purely deter¬ 
mined by its partial dissolution in the liquid ocean. 

Note that the occurrence of convection and plate tec¬ 
tonics within the ice mantle is quite uncertain. Studies by 


Eu et al. (2010) or Levi et al. (2014) suggest the possibility 


of ice mantle convection on ocean planets under certain con¬ 
ditions. Convection depends on details like the temperature 


prohle within the ice mantle, the planet’s interior heat loss 
over time, or the viscosity of the high-pressure ice phases. 
Given the uncertainties in, for example, the experimental 
data of the ice viscosity or the detailed evolution of the 
planetary interior, it is currently not possible to constrain 
the occurrence of a convective ice mantle over the lifetime of 


the planet. Additionally, it was shown by Levi et al. (2014) 


that methane (CH 4 ) could be transported through the con¬ 
vective ice mantle in the form of clathrates. To our knowl¬ 
edge, no theoretical or experimental data are available on 
clathrate formation for CO 2 under the conditions encoun¬ 
tered on an ocean planet. Assuming an efficient transport of 
carbon dioxide clathrates through the convective ice man¬ 
tle would, in principle, connect the oceanic and atmospheric 
CO 2 with the CO 2 reservoir below the ice mantle. It is a pri¬ 
ori not clear, if such a connection would work as a long-term 
stabilisation (as on Earth) or would even further enhance the 
CO 2 effect discussed in this study. Due to these uncertain¬ 
ties, we assume that no transport of CO 2 through the ice 
mantle occurs on the ocean planets in our study. 


In this study we investigate the interaction of the ocean 
with the atmospheric CO 2 . We present a potentially un¬ 
stable CO 2 cycle which, unlike the carbonate-silicate-cycle 
on Earth, acts as a destabilising climate feedback and might 
limit the overall habitability of an ocean planet. Oceanic dis¬ 
solution of CO 2 is also a well known effect in Earth climatol¬ 
ogy. The oceanic water of Earth has an overall CO 2 content 
of about two bar ( Pierrehumbert| 2010 ), which is much more 
than its actual corresponding atmospheric partial pressure 
of roughly 10~^ bar. The increased anthropogenic release 


of CO 2 into the atmosphere results in an ocean acidihca- 
tion which will considerably lower the pH value of the Earth 
ocean within the next centuries ( Caldeira Wickett||2003 ). 
While this will have a strong effect on the maritime bio¬ 
sphere in the near future, the overall impact on the atmo¬ 
sphere and surface temperature is limited to a few degree 
kelvin because of the small size of the Earth oceans (com¬ 
pared to ocean planets) and the long-term stabilisation via 
the carbonate-silicate cycle. 


While ocean planets can in principle be more massive, 
we limit our present study to planets with the mass and 
radius of the Earth, orbiting a Sun-like star, as hrst exam¬ 
ple. Otherwise, we would need additional assumptions on 
the bulk composition and internal structure of the ocean 
planet to derive a corresponding mass-radius relationship. 
This, however, is not the topic of this study which focusses 
on the impact of the temperature-dependent dissolution of 
CO 2 in the ocean on the planetary habitability. The results 
of this study are also directly applicable to ocean planets 
with higher masses and radii, though. We assume that any 
initially accumulated hydrogen and helium is already lost, 
which is a reasonable assumption for low-mass planets (AL 
ibert||2Q14 Lammer et al.||2008 ). 


The details of the model for the oceanic CO 2 dissolu¬ 
tion are described in Sect. | 2 ] In Sect. |3] we discuss the at¬ 
mospheric model and employ it to study the habitability of 
ocean planets 
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Table 1. Simon-Glatzel equation parameters for the melting 
curves of water ices. 


Ice 

Triple point parameters 

Power-law parameters 

phase 

Po (MPa) 

To (K) 

a (MPa) 

c 

Ice VT 

618.4 

272.73 

661.4 

4.69 

Ice VII2 

2216 

355 

534 

5.22 


^ [Choukroun &: Grasset| ( |2007| ^ |Alibert| ( |2014| 


2 OCEANIC CO 2 DISSOLUTION 

Both the amount of liquid ocean water (i.e. the size of the 
global ocean) as well as the dissolved CO 2 within the ocean 
(and thus also the atmospheric CO 2 content) are in general 
a function of temperature. Thus, we derive the size of the 
liquid ocean and describe the chemistry for the dissolution 
of CO 2 in ocean water in this section. We thereby do not 
account for loss processes of atmospheric CO 2 to space or 
for the potential formation of CO 2 -H 2 O clathrates as both 
of them are hard to quantify without more additional as¬ 
sumptions. 


2.1 Ocean mass 


In the following, we assume that the ocean is isothermal with 
a temperature equal to the temperature at the ocean’s sur¬ 
face to study the lower limit of the CO 2 feedback impact on 
the atmospheres of low-mass ocean planets. This assumption 
results in the smallest possible ocean size. If, for example, 
an adiabatic temperature profile had instead been adopted 
( Leger et al.||2Q04 ), the ocean would have been larger and 
able to store a larger quantity of CO 2 . This on the other 
hand would amplify the proposed CO 2 feedback mechanism 
as discussed below. 

Assuming an isothermal ocean, its maximal mass is con¬ 
strained by the maximum pressure of liquid water, given by 
the temperature and pressure-dependent melting curves of 
water ice phases VI and VII. This maximum pressure p is 
described as a function of temperature T by the empirically 
derived Simon-Glatzel equation ( Simon V Glatzel[|1929 ) 


■ Po + a ( ( ^ 


( 1 ) 


where po and To are reference values (usually the corre¬ 
sponding triple points) and a and c are power-law param¬ 
eters obtained from experimental data. The corresponding 
parameters for the equation are summarised in Tableland 
have been adopted from Alibert ( 2014| ) for ice phase VII and 
from [Ghoukroun Grasset (2007) for ice VI. 

The ocean mass is shown in Fig. as a function of 
(isothermal) temperature. It varies considerably from about 
25 to 100 Earth oceans between 273 K and 400 K. Thus, at 
low temperatures, less liquid water is available to dissolve 
the atmospheric GO 2 . 


2.2 


CO 2 dissolution model 


Following the description of Pierrehumbert (2010), we de¬ 
rive the atmospheric GO 2 content by taking into account 
the dissolution of GO 2 within the ocean. In equilibrium, the 



Figure 1. Mass of the liquid ocean and total atmospheric GO 2 
pressure as a function of isothermal ocean temperature. The at¬ 
mospheric GO 2 partial pressure (red line) is shown for a total 
GO 2 inventory of 10^^ kg and 10^® kg of N 2 . The vertical dot¬ 
ted line indicates the boundary between the formation of the ice 
phases VI and VII at the oceans bottom. 


atmospheric partial pressure of G 02 (g), Pco^{g)^ and the 
concentration of G 02 (aq) dissolved in the ocean (cco 2 (aq)) 
are linked by Henry’s law 

PCO^(g) = Kh{T) ■ Cc02(aq) (2) 

with the temperature-dependent Henry’s law constant 
Kh{T). Additional reactions occurring within the ocean 
convert the dissolved G 02 (aq) into hydrogen carbonate 
(HGO3) and carbonate (GOg”) 

C02(aq) + H20(1) ^ HC 03 -(aq) + H+(aq) , (3) 

HCOJ (aq) ^ coy (aq) + H+ (aq) . (4) 

The hydrogen ion H"*" is additionally linked to the hydroxide 
ion OH“ via the dissociation reaction of water 


H 20 ( 1 ) ^ H+(aq) + OH“(aq) . (5) 


Assuming chemical equilibrium, we solve this system 
under the constraint of charge balance 

[H+] = [HCO3-] + 2 [coy] + [OH-] (6) 


by using a Newton-Raphson method, iterating the mole 
concentration [H^] . Rate constants for the chemical reac¬ 
tions and Henry’s law have been taken from |Pierrehumbert| 
(2010). The temperature dependence of the reaction rates 
are described via the Arrhenius equation. The rate coeffi¬ 
cients used here are defined for an Earth-like oceanic salinity 
of 30%. The salinity can have a large effect on the rates of 
the reactions!^ andTests with low-salinity rate constants 


(taken from the supplementary material of Pierrehumbert 
( | 2010 | )) indicate, that our results are mostly insensitive to 
the oceanic salinity. In the cases studied in this paper, most 
of the carbon is stored as dissolved GO 2 in the ocean, such 
that the reactions and are only of minor importance. 

Note that [H^, when expressed in units of moles per 
litre, is associated with the pH value of the ocean via 


pHo 


= -logic [H" 


(7) 


Given a total GO 2 mass inventory mco^ , a second 
Newton-Raphson method is applied to find the atmospheric 
mass of GO 2 (^co 2 (g)) which the conservation of mass 

( 8 ) 
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is achieved, where mco^{a.q) is determined by the chemi¬ 
cal system described above as a function of Pco^(g)- Atmo¬ 
spheric mass and partial pressure of C 02 (g) are related via 
the condition of hydrostatic equilibrium. 

We do not consider the presence of positive calcium ions 
(Ca^^) in the charge balance. On planets like Earth, these 
ions can originate from the dissolution of limestone (CaCOg) 
in water 

CaCOg (s) ^ Ca^+ (aq) + CO^ (aq) . (9) 

Since, however, the liquid water on the planets under con¬ 
sideration is not in contact with a rocky surface, this source 
for the Ca^^ ions does not exist. Note that adding Ca^^ ions 
via dissolved limestone would strongly increase the capabil¬ 
ity of the ocean to take up CO 2 and, thereby, enhance the 
CO 2 feedback cycl^ 


2.3 Results 

Figure shows the resulting atmospheric partial pressure 
C 02 (g) as a function of isothermal ocean temperature for 
a total inventory of 10^^ kg as a particular example for an 
intermidate value of the range of CO 2 inventories considered 
in the following. The atmospheric mass of nitrogen (N 2 ) is 
fixed at an approximately Earth-like value of 10^^ kg. For 
this example, the surface pressure of C 02 (g) increases from 
about 9 bar at 273 K to more than 30 bar at 400 K. Taking 
into account the results of planetary population synthesis 


models from e.g. Marboeuf et al. (2014) to constrain the 


overall amount of CO 2 accumulated from the protoplane¬ 
tary disc for low-mass planets which have migrated into the 
habitable zone, we vary the total CO 2 inventory between 
10^® kg and 10^^ kg. The results are presented in Fig. [ 2 ] 
where a N 2 inventory of 10^® kg is used again in each case. 

We note that the results in Fig. for very high total 
CO 2 inventories are physically unrealistic as they predict 
atmospheric partial pressures of more than 100 bar. Such 
a large amount of C 02 (g) would partially condense onto 
the surface as liquid CO 2 . According to the vapour pres¬ 
sure curve ( Ambros^|1956| , the highest partial pressure of 
C02(g) at 400 K is roughly 600 bar while at 273 K only 
about 34.7 bar can be stable in the atmosphere in gaseous 
form. The condensation limit is indicated by the dotted, 
black line in Fig. Atmospheric partial pressures of CO 2 
above that line are thermodynamically unstable and would 
condense out, removing it from the atmosphere (see also 
von Paris et al. 2013a). Additionally, some of the results 
shown in Fig lie above the critical point of CO 2 , which is 
located at a temperature of 304.25 K and pressure of 73.8 
bar. Above the critical point, the description of the equi¬ 
librium chemistry from Sect. |2.2| is probably unsuitable and 
the corresponding results should be treated with great care. 
The critical point limit is indicated by the dotted, white line 
in Fig. We note that none of these unrealistic values are 
used for the atmospheric model calculations in Sect.j^ 

In general, the C02(g) partial pressure varies over sev¬ 
eral orders of magnitudes across the parameter space. Even 
though, the ocean is smaller at low temperatures and, thus. 


^ This is of course also true for all other influences affecting the 
charge balance in Eq. 



Figure 2. Atmospheric partial pressure of CO 2 as a function of 
ocean temperature and total CO 2 inventory. The amount of N 2 
in the atmosphere is flxed at 10^® kg in each case. The solid, 
black lines represent contour lines for constant CO 2 pressures. 
The condensation limit for CO 2 is indicated by the dotted line. 
Atmospheric partial pressures of CO 2 above that line are thermo¬ 
dynamically unstable and would condense out. The white, dotted 
line marks the critical point limit of CO 2 . Above this limit, the 
treatment of the equilibrium chemistry becomes unsuitable. 


less liquid water is available to dissolve the atmospheric 
CO 2 , the temperature dependence of the dissolution is 
stronger. The ability of the ocean to dissolve CO 2 decreases 
with increasing temperature, which leads to a positive, and 
thus potentially destabilizing, feedback cycle within the at¬ 
mosphere. The partial pressure of C02(g) likewise decreases 
with temperature owing to the ability of the ocean to dis¬ 
solve a larger amount of CO 2 at lower temperatures. This, 
again results in a positive feedback cycle. 

Thus, at low surface temperatures, where CO 2 would 
be required as an atmospheric greenhouse gas to heat the 
planetary surface, the ocean binds an increasing amount of 
CO 2 and, thereby, removes the important greenhouse gas 
from the atmosphere. For higher temperatures on the other 
hand, less CO 2 can be dissolved in the ocean and starts to ac¬ 
cumulate in the atmosphere, contributing to the greenhouse 
effect and further increasing the surface temperature. In con¬ 
trast to the stabilising carbonate-silicate cycle operating on 
Earth, this therefore induces positive feedback cycles, which 
could lead to potentially unstable situations: an atmospheric 
freeze-out with a runaway glaciation or a runaway green¬ 
house effect, for example. This is directly comparable to 
the water vapour feedback where the atmospheric content of 
H 2 O increases strongly with surface temperature, which may 


finally result in a runaway greenhouse process (Komabayashi 


1967 Kasting 1988). These positive CO 2 feedback cycles will 


obviously have important implications for the potential hab¬ 
itability of an ocean planet. This effect has also been noted 
in the appendix of Wordsworth & Pierrehumbert (2013) as 
a possibly important mechanism on volatile-rich planets. 

The highest surface temperatures considered in this 
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study are limited to 400 K, assuming that this is the limit 
for microorganisms to survive and thus a potential limit for 
the presence of life as we know it (|Cockell||1999| iLammer 
et al.||2009 


Rothschild & Mancinelli||2001 ). The lowest sur¬ 


face temperature for the limit of potential habitability is 
considered to be the freezing point of water (273 K). Note 
that 3D studies show that the global mean surface tempera¬ 
ture of 273 K may still allow for liquid water in the equatorial 
region (see e.g. Wolf Toon||2013 Kunze et ah poiil l 


3 ATMOSPHERIC MODEL CALCULATIONS 

Since the planetary surface temperature depends - amongst 
other things - on the atmospheric amount of CO 2 which, as 
described above, is a function of the oceans temperature, we 
use an atmospheric model coupled with the C 02 (aq) disso¬ 
lution model to constrain the potential habitability of the 
ocean planet. The one-dimensional, radiative-convective at¬ 
mospheric model was developed for terrestrial planets and 
is described in the next subsection. The results from the 
model calculations using the dissolution model and the at¬ 
mospheric model are shown in Sect. |3.2[ Focussing on the 
CO 2 feedback mechanism, we present the model calculations 
for a G2V host star: the Sun. However, this study can of 
course be applied to other central stars straightforwardly. 
The atmospheric composition for this study is restricted to 
a mixture of water, carbon dioxide, and molecular nitrogen. 
The atmospheric amount of N 2 remains fixed at an approxi¬ 
mately Earth-like value of 10^® kg. Nitrogen, unless present 
in very high amounts, has a very limited effect on the result¬ 
ing surface temperature ( von Paris et al.|2013b Wordsworth 
fc Pierrehumbert|2013 ). 


3.1 Atmospheric model description 

For the atmospheric model calculations we use a modi¬ 
fied version of the one-dimensional radiative-convective at¬ 
mospheric model extensively described in [Kitzmann et ah] 
(2010) and von Paris et al. (2010). The model has been up¬ 


dated in Stracke et al. (2015) to include new absorption coef¬ 


ficients in the shortwave part of the radiative transfer based 
on the HITEMP2010 database ( Rothman et al.|| 2010 ). The 
model considers CO 2 , H 2 O, and N 2 as atmospheric gases, ne¬ 
glecting other, potentially important, radiative trace species 
which might be present, such as O 2 , SO 2 , O 3 , or CH 4 . 

Nitrogen and CO 2 are assumed to be well mixed 
throughout the atmosphere. Water vapour profiles for low 
surface temperatures are calculated by the fixed relative hu¬ 


midity profile of Manabe & Wetherald (1967) through the 


troposphere. For the water-rich atmospheres at high stel¬ 


lar insolations, a relative humidity of unity is used (Kasting 


et al. 1993 Kopparapu et al. 2013). Above the cold trap. 


the water profile is set to an isoprofile using the H 2 O con¬ 
centration at the cold trap. 

The temperature profile is calculated from the require¬ 
ment of radiative equilibrium by a time-stepping approach. 


as well as performing a convective adjustment (Manabe & 
Strickler|1964 ), if necessary. The convective lapse rate is as¬ 


sumed to be adiabatic, taking into account the condensation 
of H 2 O and CO 2 ( Kasting et ah] 19"^ von Paris et al.|2010 ). 


As described in von Paris et al. 


radiative transfer of the atmospheric model is split into 
two wavelength regimes: a shortwave part, dealing with the 
wavelength range of the incident stellar radiation and the 
longwave part for the treatment of the atmospheric ther¬ 
mal radiation. Both parts use the correlated-/c method to 
describe the gaseous opacity. For the shortwave radiative 
transfer, 38 spectral intervals between 0.238 /xm and 4.55 /xm 
are used, in the longwave part 25 bands from 1 /xm to 500 
fim ( von Paris et al.||2010 ). In this wavelength region, only 
absorption by CO 2 and H 2 O is taken into account. Rayleigh 


scattering is considered for CO 2 , N 2 , and H 2 O (von Paris 


et al. 2010). Continuum absorption is taken into account 


for CO 2 in the longwave (von Paris et al. 2010) and H 2 O 
in both wavelength regions ( Clough et al.||19^ ). The equa¬ 
tion of radiative transfer is solved by two-stream methods 
( Toon et al.||1989 ). In the shortwave region, a (5-Eddington 
two-stream method is employed, while in the longwave part 
a hemispheric-mean two-stream radiative transfer is used 
( Kitzmann et al.|2010 ). The spectrum of the Sun describing 
the incident stellar radiation has been compiled from data 
published by [Gueymard (2004) (see Kitzmann et al. (2010) 
for details). A planetary surface albedo of 0.06 for liquid 
water is used throughout this study ( |Li et al.]|2006| . 


3.2 Results 

In principle, it is possible to couple the chemical CO 2 dis¬ 
solution model from Sect. 1 2 . 2 ] directly with the atmospheric 
model to calculate the atmospheric content of CO 2 interac¬ 
tively every time step as a function of surface temperature. 
However, this approach does not reliably converge in every 
case since one of the main atmospheric greenhouse gases 
is changed constantly within a positive feedback cycle. This 
easily leads to runaway scenarios, where CO 2 either accumu¬ 
lates in the atmosphere or is removed until the atmosphere 
freezes out. Instead of a direct coupling, we use a slightly 
different approach in this study. 

Since the desired surface temperatures are known a pri¬ 
ori (273 K and 400 K), the atmospheric content of CO 2 is 
fixed at the corresponding value from the CO 2 dissolution 
model shown in Fig. With the fixed CO 2 atmospheric 
partial pressure, the atmospheric model is used to find the 
corresponding orbital distance where the atmosphere is in 
thermal equilibrium for the chosen surface temperatures and 
CO 2 inventories. 

The resulting orbital distances for surface temperatures 
of 273 K and 400 K are shown in Fig. We additionally 
compute the distances for different stellar ages, taking into 


account the change of the solar luminosity over time (Gough 
1981). The results depicted in Fig. [^suggest that the range 
of distances for which the ocean planet remains habitable 
over an extended period of time (6 Gyr) is quite small. It is 
strongly restricted to GO 2 inventories between 10 ^® kg and 
10^^ kg and to orbital distances between 1.1 au and 1.5 au. 
Thus, without the stabilising carbonate silicate cycle, the 
habitable zone of ocean planets might be much smaller than 


predicted by the concept of the classical HZ ( [Kasting et al. 
ITopparapu et al.|2013 ). 


1993 


At large orbital distances the results are similar to the 


well-known maximum GO 2 greenhouse effect (Kasting et al. 
1993 Kopparapu et al.|2013 ). In this case, a maximum par- 


( 2010 ), the cloud-free tial pressure of GO 2 exists where the greenhouse effect still 
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Figure 3. Orbital distances for habitable surface conditions. 
Coloured lines represent surface temperatures Tg^rf of 273 K 
(blue) and 400 K (red). Line styles indicate different stellar ages of 
the Sun. The green area shows the distances and CO 2 inventories 
for which the planet remains habitable over the entire time period 
(6 Gyr). The corresponding inventories of C02(g) for Venus and 
C02(g)+C02(aq) for Earth are indicated at the right-hand axis 
for comparison. 


outweighs the Rayleigh scattering, yielding a net heating. 
For larger partial pressures, however, Rayleigh scattering 
dominates. Thus, the planet must be located closer to the 
host star to keep the surface temperature above the freezing 
point of water. At a temperature of 273 K, the total amount 
of C 02 (g) which can be present in the atmosphere is re¬ 
stricted to about 34.7 bar (equivalent to roughly 3.7 • 10^^ 
kg of total CO 2 ), because any additional CO 2 in excess of 
that partial pressure would condense onto the surface, form¬ 
ing a layer of liquid CO 2 floating on top of the ocean. 

For small CO 2 inventories, the planetary conditions at 
inner distances are mostly determined by the usual water 
vapour feedback cycle which might eventually lead to a run¬ 
away greenhouse effect (however, see Stracke et al. |2015[ for 
new results concerning the runaway greenhouse effect). In 
contrast to that, the atmospheric partial pressure of CO 2 
for larger inventories is sufficient to increase the surface 
temperature above 400 K already far beyond the current 
Earth orbit. The positive feedback cycle results here in an 
atmospheric CO 2 content larger than 30 bar, rendering the 
planetary surface uninhabitable. These findings imply that 
ocean planets are most likely continuously habitable only 
for very specific conditions. The positive feedback cycle of 
the oceanic CO 2 dissolution severely constrains the potential 
habitability of these types of planets. 

The dissolution of CO 2 in the ocean also has a strong 
impact on its pH value. Figureshows the resulting pH val¬ 
ues for the two different surface temperatures as a function 
of total CO 2 inventory. With increasing CO 2 content, more 
C 02 (aq) is dissolved in the ocean which considerably lowers 
the pH value by the release of H^ ions needed to meet the re¬ 
quirement of charge balance (Eq. ®). Starting at pH values 
of about six, it decreases to values of roughly 2.5 for the con¬ 
sidered CO 2 inventories. This may, additionally to the tem¬ 
perature effect, have an impact on the potential habitability 
of the liquid ocean ( Rothschild <V Mancineffi|2001| hammer 


et al.||20Q9 ). The pH values for surface temperatures of 400 


K are overall lower compared to the one for 273 K. This is 
related to the fact that at 273 K - even though the liquid 


Figure 4. Oceanic pH values as a function of total CO 2 inventory. 
The red line denotes surface temperatures of 400 K, the blue one 
273 K. 


ocean is smaller - more CO 2 is dissolved in the ocean. Note, 
that even at very low CO 2 inventories, the ocean of the wa- 
terworld has already a much lower pH value than the Earth 
ocean, with a pH value of about eight. This is caused by the 
absence of limestone in contact with the ocean planet’s liq¬ 
uid water which, when dissolved, would reduce the acidity 
of the ocean, buffering the pH ( Pierrehumbert||2010 ). 


4 SUMMARY 

In this study we investigated the potential habitability of 
extrasolar ocean planets. We focused on the CO 2 cycle where 
the atmospheric CO 2 can be partially dissolved in the deep 
water ocean. The high pressure at the ocean’s bottom results 
in the formation of high-pressure water ices which separate 
the liquid water from the planetary crust, thus rendering the 
usual, stabilising carbonate-silicate cycle inoperable. 

We used a chemical equilibrium chemistry to calculate 
the temperature-dependent dissolution of atmospheric CO 2 
in the liquid ocean, taking also the change of the ocean’s 
size with temperature into account. The results from the 
dissolution model imply that the CO 2 cycle operating on 
ocean planets is positive and, thus, potentially destabilis¬ 
ing. The ability of the ocean to dissolve CO 2 increases with 
decreasing temperature. Thus, at low surface temperatures 
the ocean can bind an increasing amount of CO 2 which re¬ 
moves the important greenhouse gas from the atmosphere 
where it is needed to retain habitable conditions on the sur¬ 
face. For higher temperatures on the other hand, less CO 2 
can be dissolved in the ocean and starts to accumulate in 
the atmosphere which results in the build-up of a massive 
CO 2 dominated atmosphere. This greatly contributes to the 
greenhouse effect, further increasing the surface tempera¬ 
ture. 

The results of the dissolution model were then used in 
an atmospheric model for exoplanets to calculate the dis¬ 
tances of the ocean planets to their host star where the sur¬ 
face still remains habitable for an extended period of time. 
We showed that the range of distances is very restricted in 
terms of the total CO 2 inventory of an ocean planet. 

These findings imply that ocean planets are most likely 
continuously habitable only for very specific conditions. 
The positive feedback cycle of the oceanic CO 2 dissolution 
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severely constrains the potential habitability of these types 
of planets. While this study is focused on carbon dioxide, the 
same temperature dependent solubility and therefore posi¬ 
tive feedback cycles can also be obtained for other important 
greenhouse gases, such as methane, for example. Thus, our 
results have important implications for the general question 
of habitability for exoplanets because the classical concept 
of the habitable zone for Earth-like planets which implicitly 
assumes an operating carbonate-silicate cycle is obviously 
not applicable to water-rich exoplanets. We note, however, 
that other - potentially important factors - might influence 
the results presented in this study. This includes e.g. the 
potential formation of H 2 O—CO 2 clathrates or the trans¬ 
port of CO 2 through the ice layer via clathrate diffusion or 
convection processes within the ice mantle all of which can 
affect the CO 2 content in the ocean and the atmosphere. 

Finally, we also want to point out another possible con¬ 
straint for the habitability of ocean planets. Since these plan¬ 
ets migrate from beyond the protoplanetary disks ice line, 
the water can be initially in a frozen state ( Leger et al.|2004 ) 
and has to be melted before any secondary atmosphere can 
be formed. Water ice, however, has a high albedo and, thus, 
reflects a large fraction of the incident radiation from a solar- 
type star. In the case that the primordial hydrogen-helium 
atmosphere ( Pierrehumbert Gaidos|| 2011 | , or other pro¬ 
cesses, such as impacts, the energy acquired during the ac¬ 
cretion phase within the disk, or a CO 2 /CH 4 atmosphere 
outgassed through the ice mantle ( Levi et al.|2014 ), did not 
result in an initial melting of the water ice, the stellar insola¬ 
tion must then be very high to melt the icy surface and form 
a liquid ocean. Assuming a conservative ice albedo of 0.5 and 
present day solar luminosity, the planet must be located at 
about 0.74 an to attain a global mean surface temperature 
of 273 K. At this distance, however, the planet, once the 
liquid ocean starts to form, would presently undergo a H 2 O 


runaway greenhouse scenario (Komabayashi 1967 Kasting 


1988), rendering it uninhabitable. Being a potential strong 


constraint on the habitability of ocean planets, this process 
warrants further investigation, especially with regard to the 


ice albedo feedback for different stellar types (von Paris et al. 


2013c). 
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